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Abstract 

Some observations suggest that solar spicules show small amplitude and high fre- 
quency oscillations of magneto-acoustic waves, which arise from photospheric granular 
forcing. We apply the method of MHD seismology to determine the period of kink waves. 
For this purposes, the oscillations of a magnetic cylinder embedded in a field-free envi- 
ronment is investigated. Finally, diagnostic diagrams displaying the oscillatory period 
in terms of some equilibrium parameters are provided to allow a comparison between 
theoretical results and those coming from observations. 

PACS: 96.60.-j 
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Spicules are jet-like chromospheric structures and are usually seen all around the limb of 
the Sun arising in different directions. The mechanism of spicule formation and evolution 
is not well understood (for the propulsive mechanisms, see reviews of Sterling 2000; Lorrain 
and Koutchmy 1996; Filippov et al. 2006). Spicules are relatively homogeneous in height 
along their life time of approximately 5-15 min., i.e. they are short-lived and comparable to 
the photospheric granules lifetime. They have typical up flow speeds of 20-50km/s, spicules 
diameter at chromospheric layers are of the order of 200-500km. The mean number of spicules 
per supergranule cells at height of km is approximately 40 (Pataraya et al. 1990), which 
is covered about 1 percent of the Sun's surface and they are usually concentrated between 
supergranule cells. Their temperatures and density are higher than those of the surrounding 
environment, Parenti et al. (1999) estimated n e = 10 10 cm~ 3 and T e 2 x 10 5 °K in giant 
spicules see Koutchmy and Loucif, also named macrospicules, are observed over 20,000 km 
off-limb and live 40 min. (Xia et al. 2005) in case of macro-spicules. This means that 
a magnetic field of 10G or more is needed for the low- j3 (ratio of the thermal pressure to 
the magnetic pressure) conditions (Wilhelm 2000) in case of macro-spicules. Spicule usually 
have oscillation behavior, the existence of 5 minutes oscillations in spicules have been firstly 
reported by Kulizhanishvili and Nikolsky (1978) and others including spectroscopically resolved 
observations. Recently image sequences were studied by De Pontieu et al. 2003, 2004; Xia et 
al. (2005); Ajabshirizadeh et al. (2007). These oscillations seem to be related to p-modes, but 
it is evident that if spicules are driven by p-modes, crucial details about their formation are 
still missing. Clearly, not all spicular flows are periodic, whereas most photospheric oscillations 
are. In addition, the horizontal scale for amplitude coherence of p-modes (~ 8000km) is well 
beyond the width of fibrils (De Pontieu et al. 2003). On the-other-hand, oscillations in spicule 
with shorter period have been reported by Nikolsky and Platova (1971). They found that 
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spicules oscillate along the limb with a characteristic period of about 1 min. If spicules are 
formed in thin magnetic flux tubes, then the periodic displacement of the axis observed by 
them was probably due to the propagation of kink waves. More recently, Kukhianidze et al. 
2006 have reported the observational signature of propagating kink waves in spicules. The 
period of waves was estimated to be 35-70s for a spicule with height 3500 km which may 
carry photospheric mechanical energy into the corona. The cutoff period of kink waves due to 
stratification in the hydrostatic photosphere is ~ 660s so the expected period of kink waves 
is well below the cutoff value (Singh and Dwivedi, 2007). The wavelength was found to be 
~ 1000km at the photosphere level which indicated a granular origin of the waves. Magnetic 
flux tubes support transverse kink waves that can be generated in photospheric magnetic flux 
tubes through buffeting action of granular motions (Roberts, 1979; Hasan and Kalkofen, 1999) 
although the extrapolation of the network magnetic field toward the corona is still a matter 
of discussion. In this study, the effect of gravitational stratification has been ignored, Singh 
and Dwivedi 2007 have considered kink mode periods with this effect and confirms with our 
results. 

2 Basic MHD Equations 
2.1 The Dispersion Relations 

Magneto- hydrodynamics (MHD) is one of the key tools to understand the hydrodynamics of 
magnetized plasmas. It concerns virtually all phenomena observed in the solar atmosphere: 
coronal loops, filaments, spicules, etc. Thanks to high spatial resolution, image processing, and 
time cadence capabilities of the SoHO and TRACE spacecraft, oscillating loops (and spicules) 
and propagating waves have been identified and localized in the transition region (TR) and 
chromosphere, and studied in detail since 1996. Using seeing free observations, they evidently 
complement what has been studied for a long time at ground-based, including spectroscopic 
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analysis. These discoveries established a new discipline has become know as solar atmosphere 
seismology. Many astrophysical plasmas are characterized by a set of equations that is called 
ideal MHD equations and includes the continuity the momentum, Maxwell's equations and 
Ohm's law. To understand the various oscillations and waves we observe in the spicule plasma 
we have to find wave solutions of the MHD equations (Roberts, 1981), the existence of wave 
solutions is generally derived by introducing a small perturbation of physical parameters (such 
as: density, velocity and magnetic field) of the plasma, and to derive dispersion relations u(k), 
which tell us either the group velocity or phase speed of wave. The cylindrical flux tube 
appearance of many magnetic structures in the low-/3 plasmas of the magnetosphere and more 
specifically the solar chromosphere and corona encourages an investigation of propagation in 
cylindrical geometries. Edwin and Roberts (1983) consider a uniform cylinder of magnetic 
field Bqz confined to a region of radius 2b, surrounded by a uniform magnetic field B e z;, the 
gas pressure and density within the cylinder are Pq, po , outside P e ,p e respectively (see figure 
1). The Fourier form of the velocity disturbance v± in cylindrical coordinate is: 

v\ — vi(r) exp[i(wt + n9 — k z z)\, (2.1) 

where n is integer (n=0, 1,2) which describes the azimuthally behavior of the oscillating tube 
i.e. the cylindrically symmetric (sausage or pulsation mode given by n=0, the asymmetric 
(kink or taut-wire) mode given by n=l and higher mode with n=2, 3 are called the fluting or 
interchange modes. The general governing equation is (Edwin and Roberts 1983): 

^[ " ( tli|r 2) ^ (m)] - p°wwt" - w > = °' (2 - 2) 

where m and m e are defined by m a (with a = or a = e inside or outside of the tube 
respectively): 
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where C a = (^)^ and Va<x = Ba 1 are the sound and Alfvn speed inside (or outside) the 
cylinder, and Cx a is defined as: 

C Ta = . CaVA " (2.4) 

( 7 is the ratio of specific heats.) The external and internal solutions of MHD equations need 
to be matched at the boundary by the continuity of pressure and the perpendicular component 
of velocity. After some algebra one gets the dispersion relation for magneto-acoustic waves in 
a cylindrical magnetic flux tube is found to be (Edwin and Roberts 1983; Daz et al. 2002): 

Pe (u> - fcMJ^ofgSg + P0(- 2 " ^MoWgfeg = 0, (2.5) 

where I n and K n are modified Bessel functions of order n, with I' n and K' n being the derivatives 
with respect to the argument x. This dispersion relation describes both surface (for m\ > ) 
and body waves (for ml < ). 

2.2 Kink-mode period 

Magneto-acoustic oscillations of kink mode have now been directly observed in Ha line using 
53-cm coronagraph the Abastumani Astrophysical Observatory at different heights above the 
photosphere (Kukhianidze et al. 2006). The ratio of the spicule width 2b to the spicule full 
length 2L is j = 0.01 — 0.4 , which is correspond to the dimensionless wave number (kL). 
Therefore, the observed kink- mode oscillations are in the long- wavelength regime of KL « 1, 
where the phase speed of kink-mode is practically equal to the kink speed (Spruit, 1981; 
Roberts et al. 1984) given by: 

c = , p*vX Q +p*v% e \\ (2.6) 

In the low-/? plasma limit and for the field free environment, the expression for the kink speed 
Ck reduces to: 

C k « (j^^Vao, (2.7) 
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If we denote the full spicule length 1=2L, the wavelength of the fundamental standing wave 
is the double spicule length (due to the forward and backward propagation) i.e., A = 21 and 
thus the wave number of the fundamental mode (N=l) is k z = ^ = j, while higher harmonics 
(N=2, ) would have wave numbers k z = N(j), then the time period P of a kink-mode oscillation 
at the fundamental harmonic is: 

and for higher harmonics, 

p^ ( a±^i (2.9) 



3 Results 

Such as has been pointed out before, the determinant coming from dispersion equation (2-5) 
must be truncated by taking into account a finite number of basis functions and we will use 
surface wave (m^ > in Eq. 2-5). As we know, the kink wave is essentially non- dispersive and 
has a phase speed equal to the kink speed, v p h = jr ~ Ct ~ Vao so we introduce dimensionless 
frequency which is given by ^ for odd modes and the fundamental modes and their harmonics 
have a cutoff frequency for odd modes equal ir (see Diaz et al. 2002). 

In figure 2 the eigenfrequencies of the kink (n=l) odd modes {oj cutoff — ^ ) have been plotted 
in term of the spicule half-thickness. In this plot we can see that for small value of j only the 
fundamental and lower harmonics could appear and for large ratio of -| the frequency of the 
fundamental mode is insensitive to the spicule thickness. One of the essential results of these 
plots is that the oscillatory frequency is quite insensitive to the exact value of the ratio -jr (the 
ratio of the spicule thickness to the half-length of magnetic field lines inside the spicule), i.e., 
for a given length 2L, spicules with different thickness oscillate with the same frequency, this 
is perhaps an interesting result, and this subject have been firstly reported by Daz et al. 2002, 
for kink waves in the prominence fibrils. We continue by choosing two typical values of the 
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spicule full length, namely, 1=3500 and 8500km, and then plot the frequency of odd modes 
below the cutoff as a function of (^) (see figure 3-a, b). 

In these plots, we keep j = 0.05, 0.03 constant for two lengths of spicule respectively. For 
all harmonic modes the frequency is seen to slowly decrease with increasing (^) and more 
oscillatory modes are present for large values of this quantity. To compute the corresponding 
periods for each harmonic modes, we used 1=3500 and 8500km, B = 30G and p = 3 x 
I0~ w kgm~ 3 . The periods obtained from these quantities and using Eq. (2-9) is labeled P and 
correspond to the right vertical axes in figure 3 (see also Table 1). 

The period, P, (figure 3- a, b right vertical axes and table 1) for different values of \ , L 
calculated from the expression of Eq. (2-5) and (2-9), show that the fundamental and first 
harmonic have typical periods of ~ 30 — 80s. These periods which are obtained from a simple 
MHD approach are found in the more recently reported observational results (Kukhianidze et 
al. 2006). where oscillation periods in the range of kink wave of ~ 35 — 70s are found for 
spicules with height 3500-8700 km. From the theoretical point of view, we therefore expect 
kink wave periods in the range of 80 — 120s for length of ps 10, 000 — 14, 000km which could 
be found from future observations of kink wave inside spicules (Hinode observations). 
Acknowledgments. The authors are most grateful to A. J. Daz for useful discussions and 
critical reading of the manuscript. 
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Figure Captions 

Figure-1: Fig.l- Sketch of the equilibrium configuration used in this study. The density 
and magnetic field inside the spicule are B e , p e , and in chromosphere environment are B ,p . 
The magnetic field is uniform and along z-axis.. 

Figure-2: Fig. 2. Variation of the frequency with the spicule half-thickness for Kink- 
Modes for the set of value 2L=8500km, and for the density ratio ^=0.03. 

Figure-3: Fig. 3- a, b. Frequency of the kink odd modes vs.^ for two full length of 
sipules and -| = 0.05, 0.03 for the parameters . The right-hand axis provides the period P after 
estimation that magnetic field strength, the spicule density, and full-length of spicules are 
p = 3 x I0~ 10 kgm~ 3 (corresponding to a 2 x 10 11 number density), B = 30G, and 2L=3500, 
8500km, respectively. 

Figure-4: TABLE 1 



